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Abstract. We investigate the microscopic origin of the relaxation time coefficient 
in relativistic fluid dynamics. We show that the extraction of the shear viscosity 
relaxation time via the gradient expansion is ambiguous and in general fails to 
give the correct result. The correct value for the shear viscosity relaxation time 
is extracted from the slowest non-hydrodynamic pole of the corresponding retarded 
Green's function, if such a pole is purely imaginary. According to the AdS/CFT 
correspondence, in strongly-coupled Af = 4 SYM the non-hydrodynamic poles of the 
shear stress tensor nearest to the origin have a nonzero real part, which implies that 
the transient fluid-dynamical equations for this gauge theory are not equivalent to the 
well-known Israel-Stewart equations. 

The derivation of a consistent theory of relativistic fluid dynamics has been a 
challenge for some time. The difficulty resides in the parabolic nature of Navier-Stokes 
theory, which allows signals that propagate with infinite speed. While in nonrelativistic 
theories this feature can be overlooked, in relativistic systems causality is intimately 
related to stability [1] and this feature leads to unstable equations of motion. A causal 
and stable theory of fluid dynamics can be obtained by taking into account the transient 
dynamics of the dissipative currents [2], i.e., by considering the characteristic times 
within which fluid-dynamical dissipative currents, such as the shear stress tensor tt^ , 
relax towards their asymptotic Navier-Stokes solution. 

In this contribution, we investigate the physical origin of such characteristic time 
scales for transient fluid behavior, called the relaxation time coefficients [2]. We 
discuss two distinct formulations of relativistic fluid dynamics: The first, derived in 
Ref. |3], extends Navier-Stokes theory via the gradient expansion, for the specific case 
of conformal fluids, and introduces the relaxation time as a type of regulator to obtain 
stable fluid-dynamical equations. This procedure is implemented in such a way that 
the regularized theory is equivalent to the gradient expansion for low frequencies and 
large wavelengths. The second approach, originally presented in Ref. [I], derives fluid 
dynamics from the underlying microscopic theory taking into account the slowest non- 
hydrodynamic mode for each dissipative current. The theory obtained in this case is 
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not only equivalent to the gradient expansion for low frequencies and large wavelengths, 
it is also able to describe the transient dynamics of the dissipative currents. In this 
contribution, we consider only the case of the shear stress tensor. 

In Ref . [3] , a relativistic Burnett equation is derived for conformal fluids. The basic 
idea of this derivation is to extend Navier-Stokes theory by including all possible terms 
of second order in gradients of temperature T and four-velocity m m that are symmetric, 
transverse, traceless, and that transform homogeneously under Weyl transformations 

®, 

Of = RW - 2 (V<" V v) InT - V <M In TV") lnT) , 

Of = - 2u a R a ^up, 

Of = <7<V> A , Of = <7<fV> A , Of = u^ x u v >\ 

with being the shear tensor, the vorticity tensor, R^ v the Ricci tensor, 
and R^ ua ^ the Riemann tensor. We adopt the notation = A^A al3 with 

A^ = (^A£Ap + A^A^ — 2/3A MI/ A Q; 3^ /2 being the traceless and symmetric projection 
operator and A^ v = g^ v + u^u v being the projection operator orthogonal to w M . Then, 
the shear stress tensor, up to second order in gradients, is assumed to be 

5 

tt"" = -2^ + ^AiOf . (1) 

i=l 

As is well-known, the relativistic and nonrelativistic Burnett equations are unstable 
and have no practical use. In order to render the theory stable, a relaxation time 
coefficient for the shear stress tensor must be introduced. In Ref. [3] this was done in 
the following way. Using the conservation laws of energy and momentum it is possible 
to show that, up to third order in gradients, 

& (^) + i a ^0 = of _ of _ iof + 20f + O (V 3 ) . (2) 
This equation is then used to replace the term Of in Eq. (1), i.e., 
tt^ = - 2 V a^ + A x + g^fl) + (A 2 + Ai) Of 

+ (A 3 + i A x ) Of + A 4 Of + (A 5 - 2Ax) Of. (3) 

The last step of the derivation is to construct a relaxation equation for tt^ by 
substituting the first-order solution, a^ v = — {l/2ri) 7rf in all second-order terms 

tvtt^ +n^= - 2 V a^ - \t^ v Q + (A 2 + Ax) Of 

+ (a 3 + Of + A 4 Of + (Ab + 2Ax) Of, (4) 

where Of corresponds to Of with the substitution a^ u — > —ix^ v /(2rj) and the 
coefficient was identified as Ai/ (2rj). 

The coefficients Ai t 2 can be calculated using linear response theory via metric 
perturbations. In strongly-coupled Af = 4 SYM theory, the coefficients Ai and r\ 
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were obtained from derivatives of the retarded Green's function at the origirJ|| and 
the relaxation time was computed using this prescription, i.e., r n = XxJ (2rj) = 
(2-1h2)/(2ttT). 

Note, however, that this prescription to obtain the shear viscosity relaxation time is 
ambiguous and, consequently, all the other transport coefficients of Eq. (Hj) are also ill- 
defined. To the best of the authors' knowledge, the substitution of the second-order term 
Oi" using Eq. fl2]) is a choice. For example, if instead of Of", the term is replaced 
using Eq. (J2J), the relaxation time would be identified as r n = A2/ (2??). Naturally, 
there are infinitely many ways to make such a replacement and, in this sense, infinitely 
different values for the shear viscosity relaxation time could have been obtained in this 
approach. We remark that any choice would satisfy the conditions set by the authors 
of Ref. [3] that the regularized theory in Eq. (j4]) should be equivalent to the gradient 
expansion at low frequencies and large wavelengths. 

The reason for this ambiguity can be easily understood. Different choices of 
replacements would describe different transient dynamics, described by their respective 
relaxation time coefficients, but would invariably lead to the same asymptotic solution, 
i.e., the Burnett equation in Eq. ([[]). Naturally, it is always possible to determine the 
asymptotic dynamics of a system once the transient dynamics is known, but the inverse 
is impossible. Since the starting point of Ref. [3] is the asymptotic theory itself, i.e., 
the gradient expansion. The actual transient dynamics with the correct relaxation time 
cannot be uniquely extracted from this approach. 

Recently, a new method to derive fluid dynamics was introduced in Ref. jl] . In this 
approach, transient fluid dynamics is derived directly from an underlying microscopic 
theory. Furthermore, this work derives a condition for the retarded Green's function 
under which the linearized equation of motion of a dissipative current can be reduced to 
a relaxation-type equation. The equations of motion obtained were shown to have the 
gradient expansion as its asymptotic solution, i.e., they are equivalent to the Navier- 
Stokes and Burnett theories at low frequency and large wavelength. However, these 
equations of motion are not only capable to describe the asymptotic regime but they 
also correctly describe the dynamics of the shear stress tensor at time scales of the 
order of the relaxation time, which was shown to be related to the pole of the retarded 
Green's function closest to the origin. This clears up a few misconceptions in the field 
such as whether the relaxation time is of microscopic origin or of fluid-dynamic origin 
and whether a truncated Taylor expansion around the origin can be used to correctly 
obtain the relaxation time of dissipative fluid dynamics. 

This approach was applied in Ref. [4] to study dilute gases where the Boltzmann 
equation can be used to describe the underlying microscopic theory. In this work, 
microscopic formulas for the shear viscosity coefficient, 77, and relaxation time, r n , were 
derived using linear response theory. It was demonstrated that the retarded Green's 
function associated with the shear stress tensor has infinitely many simple poles along 

I The associated retarded Green's function in this case was computed using the AdS/CFT 
correspondence [5]. 
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the imaginary axis making it possible to reduce the equation of motion for the shear 
stress tensor to a relaxation-type equation. The formulas for rj and r n were shown to be 
consistent with those previously derived by Israel and Stewart (2j [6] via the 14-moment 
approximation. This puts relativistic viscous hydrodynamics on a solid ground. 

For strongly-coupled conformal plasmas defined via the AdS/CFT correspondence 
the situation was shown to be more subtle [TJ. For a relaxation equation to exist, i.e., 
an equation of motion with only one time derivative, it is mandatory for the pole of 
the retarded Green's function nearest to the origin to be purely imaginary. In this case, 
the system obeys a simple relaxation equation at long times, as it is the case with the 
Boltzmann equation pE]. This, however, does not occur in strongly coupled M = 4 
SYM theory: in this theory the poles of the retarded Green's function nearest to the 
origin have real parts and are symmetric around the imaginary axis [5]. Thus, both 
poles contribute equally to the long-time dynamics. Since the poles have a real part, 
the system not only relaxes to its Navier-Stokes solution, but oscillates around it as 
well. Then, the linearized equation of motion satisfied by the shear stress tensor, at 
long times, should be of the form: 

X2T^ V) + Xi^ u) + ^ = + • • ■ • (5) 

Note that this equation was not assumed: It was directly derived from the retarded 
Green's function of the microscopic theory and it describes the long-time dynamics of 
the system. It should also be remarked that this equation has the Navier-Stokes and 
Burnett theories as its asymptotic solution. The coefficients of the equation above in 
the case of strongly coupled M = 4 SYM theory were computed in Ref. [7J and shown 
to be 

Xi = 1.27/(4ttT), X 2 = 0.93/(4ttT) 2 , V /s = 1/ (4tt) , (6) 

where T is the temperature and s is the entropy density. 
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